Anisotropic superfluidity in a dipolar Bose gas 
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We study the superfluid character of a dipolar Bose-Einstein condensate (DBEC) in a quasi-two 
dimensional (q2D) geometry. In particular, we allow for the dipole polarization to have some non- 
zero projection into the plane of the condensate so that the effective interaction is anisotropic in 
this plane, yielding an anisotropic dispersion for propagation of quasiparticles. By performing direct 
numerical simulations of a probe moving through the DBEC, we observe the sudden onset of drag or 
creation of vortex-antivortex pairs at critical velocities that depend strongly on the direction of the 
probe's motion. This anisotropy emerges because of the anisotropic manifestation of a roton-like 
mode in the system. 



A quintessential feature of a superfluid is its ability to 
support dissipationless flow, for example, when an ob- 
ject moves through a superfluid and experiences no drag 
force. This, however, only occurs when the object is mov- 
ing below a certain critical velocity; when it exceeds this 
critical velocity it dissipates energy into excitations of 
the superfluid, resulting in a net drag force on the object 
and the breakdown of superfluid flow. 

In many supcrfluids, such as dilute Bose-Einstein con- 
densates (BECs) of atoms, this critical velocity is simply 
the speed of sound in the system, which is set by the 
density and the s-wave scattering length of the atoms. 
However, for a dense superfluid such as liquid 4 He, this 
is not the case. In 4 He, the critical velocity is set by a 
roton mode, corresponding to a peak in the static struc- 
ture factor of the system at some finite, non-zero momen- 
tum, with a characteristic velocity that is considerably 
less than the speed of sound in the liquid. This fea- 
ture has been verified experimentally via measurements 
of ion-drift velocity in the fluid [l[ , thereby providing in- 
sight into the detailed structure of the system. 

Interestingly, a BEC of dipolar constituents (DBEC) is 
also expected to possess a roton-like dispersion, in spite 
of existing in a dilute gaseous state ■ Unlike the disper- 
sion of 4 He, the dispersion of a DBEC is highly tunable 
as a function of the condensate density and dipole-dipolc 
interaction (ddi) strength. Additionally, the DBEC is 
set apart from liquid 4 He in that its interactions de- 
pend on how the dipoles are oriented in space. Thus, the 
DBEC provides an ideal system to study the effects that 
anisotropics have on the bulk properties of a superfluid, 
such as its critical velocity. Anisotropic dispersions have 
been predicted for a ID lattice system of q2D DBECs Q, 
periodically dressed BECs Q and for dipolar gases in 
a 2D lattice Q. Additionally, anisotropic solitons have 
been predicted for dipolar gases [f|. 

In this Letter we consider a DBEC in a quasi-two- 
dimensional (q2D) geometry and allow for the dipoles to 
be polarized at a nonzero angle into this plane so that the 
in-plane interaction is anisotropic. We perform numer- 
ical simulations of a probe moving through the DBEC. 
This probe experiences a sudden onset of drag at a cer- 
tain velocity, the critical velocity v c , which is very dif- 



ferent depending on the direction of the probe's motion. 
For perturbative probes the drag originates from the pro- 
duction of quasiparticles and v c is nearly equal to Vl, the 
critical velocity predicted by the dispersion of the DBEC. 
In contrast, for a larger probe vortex-antivortex pairs are 
formed, usually leading to a critical velocity smaller than 
vl @, 13 • For both cases we find that the critical veloc- 
ity is larger in the direction of the dipole tilt than in 
the perpendicular direction. Interestingly, while the ro- 
ton displays an anisotropic character, the speed of sound 
in the system remains isotropic. Thus, we characterize 
the DBEC as an anisotropic superfluid while illuminat- 
ing the crucial role that the roton plays in generating this 
anisotropic behavior. Importantly, such a system offers 
a path to study a bosonic system which has anisotropic 
collective behavior which is typically found in composite 
fermion systems, such as high T c d-wave Superconduc- 
tors. 

We work in the q2D geometry by assuming that there 
is a strong one-dimensional (ID) harmonic trap in the 
z-direction, Uid(z) = ^mui z z 2 , where m is the bosonic 
mass and u z is the trapping frequency. This allows the 
condensate wavefunction to be written in the separa- 
ble form ^(r, t) = x(z)ip(p,t) where x( z ) 1S the ground 
state harmonic oscillator wavefunction. By inserting this 
ansatz into the Gross-Pitaevskii equation (CPE), de- 
scribing the dilute, zero-temperature DBEC, and inte- 
grating out the z-dependence, we derive the modified, 
time-dependent GPE for the q2D system, 



ihd ^ = { -— V + V p + 9\4'\ + 9d$ } 1>, (1) 
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where tp = tp(p,t) is the in-plane condensate wave- 
function, V p — Vp(p, t) is a time-dependent probe po- 



is the mean-field coupling 
is the scattering length, 



tential, g = 2\2~kTi a s /ml z 
for contact interactions, a s 
l z = ^Jh/muj z is the axial harmonic oscillator length, 
9d = d 2 / \/2~nl z is the ddi coupling and gd<& = gd&(p, t) is 
the mean-field potential due to the ddi, where $ (p, t) = 
-^-J : ~ 1 [h(\i, t)F(^=)}, by the convolution theorem. Here, 
J- is the 2D Fourier transform operator and n(k, t) = 
T[n(p,t)}. 
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The function ^j-gdF(q), where q = kl z /y/2, is the k- 
space ddi for the q2D geometry. It has two contribu- 
tions coming from polarization perpendicular or parallel 
to the direction of the dipole tilt, F(q) = cos 2 (a)i 7 l(q) + 
sin 2 (a)F|| (q) where a is the angle between z and the 
polarization vector d. These contributions are -Fy (q) = 
— 1 + 3y/n(qj/q)e q2 eric(q), where q^ is the wave vector 
along the direction of the projection of d onto the x- 
y plane, erfc is the complementary error function and 
Fj_(q) = 2 - 3^e? 2 erfc( g ). 

To simplify this problem, we rescale energies in units 
of the chemical potential, given by u* — on {l + 
(3^-(3 cos(a) 2 — 1)} for the unperturbed system jjj]. Here, 
f3 = gd/g- This leads to characteristic units of length 
given by the coherence length £* = h/y/mu*; time 
t* = %/fJj*] and velocity c* = y/ fi* /to. Additionally, 
we rescale the wavefunction ip — > ip/ y/n~Q where no is the 
integrated 2D density of the unperturbed system. The 
rescaled coupling constants are then g* = gn /fi* and 
9*i = Sd n o/u* = fig* ■ This formalism describes the q2D 
DBEC, including the ground (condensed) state and the 
dispersion relation that describes the system's quasipar- 
ticle spectrum. 

The dispersion relation of a homogeneous q2D DBEC 
is given in Bogoliubov theory by fioj 
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;(k) = 
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For a — (polarization along the trap axis) this disper- 
sion does not depend on the direction of the quasiparticle 
propagation. However, for a ^ 0, or for nonzero projec- 
tion of d onto the x-y plane, the direction of k becomes 
important in describing the quasiparticles of the system. 
Landau famously used the dispersion relation of a Bose 
fluid to determine its superfluid critical velocity, defining 
the Landau critical velocity, vl = min[w(k)/fc] As 
vl is given in terms of the dispersion relation w(k), it 
then also depends on the direction of k, and thus is an 
anisotropic quantity when | cos (a)\ < 1. 

To illustrate this point, we use the parameters g* = 
0.25, g* d = 0.36, l z /C = 0.5 and a = tt/4 (fi* = 1), 
which are chosen to best illustrate anisotropic effects 
while keeping safely away from the unstable regime; we 
will identify them with experimental parameters below. 
Figure [TTa) shows the dispersion calculated using these 
parameters for quasiparticle propagation parallel to (||) 
and perpendicular to (_L) the tilt of the dipoles into the 
plane. For parallel propagation, the dispersion resembles 
that of a system with contact interactions; the curve goes 
smoothly from the linear phonon regime at small k to the 
free-particle regime at large k. For this case, vl/c* = 1.0 
(dashed red line), meaning that the critical velocity is 
identical to the speed of sound. In contrast, the perpen- 
dicular dispersion curve exhibits a roton-like feature at 
intermediate k, setting vl/c* = 0.50 (dashed blue line). 
The inset in figure HJa) shows vl as a function of the az- 
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FIG. 1: (color online) (a) The Bogoliubov dispersions for g* = 
0.25, — 0.36, lz/£* = 0.5 and a = tt/4 for propagation 
perpendicular (_L) to and parallel (||) to the tilt of the dipoles, 
shown by the blue and red lines, respectively. The dashed 
lines have slopes that are the Landau critical velocities (vl) of 
the dispersions, while the inset shows vl as a smooth function 
of the angle g between the || (g/n = 0) and _L (g/n — 1/2) 
propagation directions, (b) A contour plot of the density for a 
stationary obstacle with amplitude A p /g* = 1.0. The shaded 
region indicates a density exceeding 1.05no, and the arrow 
indicates the direction of polarization, (c) Density slices of 
(b) along the parallel (dashed red) and perpendicular (blue) 
directions. The density oscillation due to the roton is clear in 
the perpendicular case. 



imuthal angle, to the angle between k^ and k. Interest- 
ingly, the speed of sound, given by c s = lim/ c _ i .o[a'(fc) / 'k], 
is the same for both parallel and perpendicular propaga- 
tion, c s /c* = 1, and is in fact isotropic. Therefore, the 
anisotropy in the spectrum occurs only at finite k due to 
the presence of an anisotropic roton. 

The impact of the anisotropic roton can be directly 
seen in the density of the gas. In Figure [1Tb) we 
show a contour plot of the density in the presence of 
a repulsive Gaussian potential, or "probe," of the form 
V p (x, y) = A p cxp (—(x 2 + y 2 )/<jp) with A p //j,* = 1.0 and 
Up/S,* = 2.0, as may be realized by shining a blue-detuned 
laser on the system from the z-direction. The shaded re- 
gions indicate density above 1.05rio, and the arrow indi- 
cates the direction of polarization. In FigureQTc) we plot 
density slices of this distribution to more clearly show the 
density profile in the parallel (dashed red) and perpen- 
dicular (solid blue) directions. Interestingly, the high- 
density regions occur in the direction perpendicular to 
the tilt of the dipoles, the same direction that exhibits a 
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roton feature in the dispersion. Indeed, it was shown in 
ref. [HI that a DBEC in the presence of a repulsive Gaus- 
sian (or a vortex core), at sufficiently large density, will 
exhibit density oscillations due to the manifestation of 
the roton. Here, we see a manifestation of the anisotropic 
roton in the static structure of the q2D DBEC. 

We now address the question of what happens to this 
anisotropic DBEC when the probe is moved through it 
with velocity v, by numerically solving Eq. ((TJ) with the 
Gaussian potential V p (x — vt,y). For concreteness, we 
consider motion parallel and perpendicular to the tilt of 
the dipoles by tilting d into the x and y directions, re- 
spectively, while fixing the direction of the probe velocity 
so that v = vx. 

Figure [21a) shows the time-averaged drag force (aver- 
aged up to t = IOOt*) acting on a "weak" probe with 
parameters A p //i* = 0.1 and c P /£* = 2.0. The force at 
time t is given by F(t) = - J d 2 p\%j;{p,t)\ 2 VV p {p,t) Q. 
In this case, the probe is sufficiently weak so that no vor- 
tices are nucleated in the fluid, and instead only quasi- 
particles are excited. The presence of a force on the probe 
signifies the excitation of quasiparticles, and thus the 
breakdown of superffiiid flow. There is a clear anisotropic 
onset of force in these simulations that agrees very well 
with the anisotropic vl given by the Bogoliubov dis- 
persions in figure [lja), resulting in critical velocities of 
v c /c* = 0.90(0.46) for parallel (perpendicular) motion of 
this probe, determined by the velocity at which the drag 
force suddenly rises. 

It has been shown that vl is recovered as the true crit- 
ical velocity only when the superfluid is perturbed by 
a vanishingly small object (13, Il4j]. Additionally, while 
quasiparticle excitations are a natural feature to study 
when considering the breakdown of superfluid flow, they 
may be difficult to observe experimentally, especially in 
the limit where the probe is perturbative. Vortices, on 
the other hand, are superfluid excitations in the form of 
topological defects that create regions of zero density and 
are easier to observe experimentally than quasiparticles. 
The first measurements of v c in a BEC were from obser- 
vations of the sudden onset of heating [l^, [l6[ , believed 
to be related to vortex production in the BECs. More re- 
cently, Ref. [13] used experimental finesse to controllably 
create vortex pairs to observe v c . 

Motivated by these circumstances, we investigate the 
critical velocity for vortex formation in the q2D DBEC by 
using a moving probe with an amplitude that is linearly 
ramped from A p = to A p //jl* = 1.0 in a time 10r* with 
<7p/£* =2.0. The critical velocity in this case corresponds 
to the probe velocity above which vortices are formed, 
signaling the breakdown of superfluidity as energy is used 
to create these topological excitations. 

We observe a significant difference in the critical veloc- 
ity at which vortices are formed between a probe moving 
parallel and perpendicular to the dipole polarization. In 
Figure [2{b) we show the maximum number of vortices 
formed within t/r* = 100. The critical velocities are 
v c = 0.46(0.28) for motion parallel (perpendicular) to 
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FIG. 2: (color online) (a) The mean drag force acting on the 
probe with A p /fj,* = 0.1 and cr p /£* = 2.0 calculated up to 
time t/r* = 100 for motion perpendicular to (blue diamonds) 
and parallel to (red squares) the dipole tilt. The dotted lines 
represent the corresponding vl, in excellent agreement with 
the numerical simulations, (b) The maximum vortex num- 
ber produced by a probe with A v j 'fx* = 1.0 and <t p /C* = 2.0 
calculated up to t/r* = 100 for motion perpendicular to (blue 
diamonds) and parallel to (red squares) the dipole tilt. The 
corresponding critical velocities are: vi^/c* = 0.27 and 
v ( J l) /c* = 0.46. 



the dipole tilt. These values are about half the value of 
the critical velocities obtained using the weaker probe, 
but this is not unexpected @,@]- I n a superfluid, vortices 
have quantized circulation: <f> v • dl = 2irhn/m, where v 
is the velocity field of the fluid and n is an integer, cor- 
responding to phase winding of 2-7rn around the vortex 
core. We count vortices in our simulations by finding the 
pha se winding on a plaquette of neighboring grid points 

The physical mechanism that sets the critical velocity 
for vortex formation is not rigorously understood. How- 
ever, it is theorized that the maximum local fluid veloc- 
ity about an obstacle, being larger than the background 
flow velocity, sets the critical velocity via the Landau 
criterion. This idea has been fruitful 0, Q, and we find 
qualitative agreement with this theory here, as the di- 
rection with lower v c is also the direction of flow most 
likely to spawn vortices. However, we note that the ddi 
is anisotropic although the fully condensed (ground) state 
of the system is completely isotropic. The anisotropics 
only appear in the dispersion relation and in the ground 
state of the system in the presence of a perturbing po- 
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FIG. 3: (color online) (a) and (b) The densities are shown for 
the || (top, red) and _L (bottom, blue) cases for v/c* =0.30 and 
t/r* = 125. Only the _L case has exceeded its v c . (c) and (d) 
Densities for v/c* =0.50 at t/r* = 100. The parallel case (top, 
red) has exceeded its critical velocity and the perpendicular 
case (bottom, blue) has been wildly excited. The obstacle is 
at x = y = and moving in the +x direction. The shaded 
regions of this plot occur when the density exceeds 1.05no, 
and the arrow indicates the direction of dipole tilt. 

tential, which is intimately related to the dispersion re- 
lation [12j. Thus, the anisotropies in the critical velocity 
for vortex formation are due to the anisotropy of the ro- 
ton mode, just like the critical velocities for quasiparticle 
excitations. 

Figure [3] shows contour plots of the condensate density 
for both parallel (red contours, top row) and perpendic- 
ular (blue contours, bottom row) motion of the probe 
relative to the dipole tilt for velocities v/c* = 0.3 (left 
column) and v/c* = 0.5 (right column), where the probe 
is moving in the x direction and is located at the origin 
at the time shown in figure Recall that v/c* = 0.3 
is just above v c for vortex formation for perpendicular 
motion, but well beneath v c for parallel motion. This is 
reflected in the figure, where in (a) no vortices have been 
formed for parallel motion, while in (b) a vortex pair has 
been formed for perpendicular motion for the same probe 
velocity. 

For the case of v/c* = 0.5, we see that the parallel case 
in (c) has formed a vortex pair, and in (d) the perpendic- 



ular case has been wildly excited. There is an important 
contrast to be made in the density profiles when there is 
a single vortex pair in (b) and (c). In the parallel case 
(c) we see that a high density region occurs between the 
vortex pair and is elongated in the polarization direc- 
tion. In contrast, for the perpendicular case (b) there is 
a low density region between the vortex pair and high 
density regions on either side of the vortex pair. Both 
the anisotropic supcrfluid critical velocity for vortex pair 
production and these contrasting density profiles present 
means to observe the effects of the roton in DBEC di- 
rectly. 

In addition to investigating this q2D system, we have 
performed simulations for a fully trapped DBEC. For a 
trap aspect ratio of A = lu z /lo p = 50, we find critical 
velocities for vortex production are strongly anisotropic, 
and the v c are numerically similar to the free case. In 
these simulations, we start the probe in the center of 
the trap and move outwards in the parallel or perpen- 
dicular direction, linearly ramping the amplitude of the 
laser down to zero by the time it reaches the zero den- 
sity region. Such a simulation is experimentally realiz- 
able in a DBEC of atomic 52 Cr, for example, having a 
permenant magnetic dipole moment of 6/is where (1b is 
the Bohr magneton, for a DBEC with particle number 
N ~ f8.5 x 10 3 , scattering length a s = 5.0ao where ao is 
the Bohr radius, radial trap frequency w p = 27r x 20Hz 
and a blue-detuned laser with width a p = 1 .76/^m. The 
speed of sound in this system is c s = 0.f6cm/s in the 
center of the trap [20| . 

In conclusion, we have characterized the DBEC as an 
anisotropic superfluid by performing numerical simula- 
tions of a blue-detuned laser moving through the system 
in directions parallel and perpendicular to the dipole po- 
larization. We find a sudden onset of drag on the laser 
at velocities that depend strongly on the direction of 
motion, and attribute the anisotropy in critical veloc- 
ity to the anisotropic roton so that a measurement of an 
anisotropic critical velocity in a DBEC corresponds to a 
measurement of the roton in the system. Additionally, 
by considering a DBEC that is experimentally realizable 
with atomic 52 Cr, we propose a single, stable constituent 
with which to study anisotropic superfluidity, while other 
systems such as supcrfluids of <i-wave Cooper pairs are 
more conceptually and experimentally difficult to control. 
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